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Abstract 



We extend a previous improved action study of the Landau gauge gluon prop- 
agator, by using a variety of lattices with spacings from a = 0.17 to 0.41 fm, 
to more fully explore finite volume and discretization effects. We also extend 
a previously used technique for minimizing lattice artifacts, the appropri- 
ate choice of momentum variable or "kinematic correction", by considering 
it more generally as a "tree- level correction". We demonstrate that by using 
tree-level correction, determined by the tree-level behavior of the action being 
considered, it is possible to obtain scaling behavior over a very wide range of 
momenta and lattice spacings. This makes it possible to explore the infinite 
volume and continuum limits of the Landau-gauge gluon propagator. 
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I. INTRODUCTION 



There has long been interest in the infrared behavior of the gluon propagator as a probe 
into the mechanism of confinement and lattice studies focusing on its ultraviolet behav- 
ior have been used to calculate the running QCD coupling 0]. In this report we use the 
propagator as a test-bed for an improved action and also as a means to investigate a general 
tree-level correction technique. 

The infrared part of any lattice calculation may be affected by the finite volume of the 
lattice. Larger volumes mean either more lattice points (with increased computational cost) 
or coarser lattices (with corresponding discretization errors). Improved actions have been 
shown to be effective at reducing discretization errors at a given lattice spacing in studies of 
the static quark potential 0] and the hadron spectrum 0] and have become a necessary part 
of finite temperature studies |^ . The desire for larger physical volumes thus provides strong 
motivation for using improved actions. We study the gluon propagator, in Landau gauge, in 
quenched QCD (pure SU(3) Yang- Mills), using the mean-field (tadpole) improved version 
of the tree-level, 0{a^) Symanzik improved gauge action f^HH- 

To assess the effects of finite lattice spacing, we calculate the propagator on a set of 
lattices from 8^ x 16 at = 3.75 having a = 0.413 fm to 16^ x 32 at /3 = 4.38 having 
a = 0.167 fm. To assist us in observing possible finite volume effects, we add to this set 
a 16'^ X 32 lattice at /3 = 3.92 with a = 0.353, which has the very large physical size of 
5.65^ X 11.30 fm^. Some preliminary results of this work were reported in Ref. . 

We will show that tree-level correction reduces rotational symmetry breaking and dra- 
matically improves the ultraviolet behavior of the propagator and thus the approach to the 
continuum limit. For lattices as coarse as 0.17 fm the gluon propagator has surprisingly 
good behavior for the entire range of available momenta. The infrared behavior of the gluon 
propagator is robust even with an extremely coarse lattice spacing of 0.41 fm. Our calcu- 
lations on a lattice with a large volume indicates that finite volume effects are small. The 
Landau gauge gluon propagator is again found to be infrared finite, in agreement with ear- 
lier studies. The combination of an improved action with appropriate tree-level correction 
appears to be a powerful tool. The generalization of these methods to the study of other 



Green's functions will be discussed in a forthcoming work [pJ]] . 

II. THE LANDAU GAUGE GLUON PROPAGATOR 

We employ the tree-level, mean-field improved gauge action of Liischer and Weisz [|,^ 

= S„.. + 0(0") + O(a^g^), (2.1) 
where P^,^ and R^^i, are the usual plaquette and rectangle operators 

Pm-(x) = U^{x)U,{x + fi)Ul{x + y)Ul{x) (2.2) 

and 
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Ri^uix) = U^{x)U^{x + fi)U^{x + z> + (i)Ul{x + 2v)Ul{x + v)Ul{x) 

+U^{x)U^{x + fi)U,{x + 2fi)Ul{x + /i + u)Ul{x + z>)t/t(a;), (2.3) 

and Nc = "i is the number of colors. We use the plaquette definition for the tadpole factor 

«o = ^l^7^eTr(P^,)V . (2.4) 



Our gauge field configurations were generated using the Cabbibo-Marinari |jl2| pseudo- 



heatbath algorithm with appropriate link partitioning |13 



Given that the gauge links U^{x) are expressed in terms of the continuum gluon fields as 

U^{x) = pe*9/o (2.5) 
the dimensionless lattice gluon field A^{x) may be obtained from 

A,{X + = ^m^) - f^;(^)}t.aceless (2-6) 

which is accurate to (!?(a^). This is, of course, only one of many possible ways to calculate 
the gluon field on the lattice. In Eq. (|2.6| ), is calculated at the midpoint of the link 
to remove 0{a) terms. Note that we have also included the tadpole factor to improve the 
normalization. 

We calculate the gluon propagator in coordinate space 

D;l{x,y)^{A';,{x)Al{y)), (2.7) 

using Eq. ( p.6|) . To improve statistics, we use translational invariance and calculate 

K(y) = T7(E^»Kix + y))- (2.8) 

The quantity that will be of interest to us is the scalar part of the propagator in momentum 
space, so first we take the trace over color components 

then sum over the Lorentz components^ of the Fourier transform 

m = ^T^j:Y.e'''D,,{y) (2.10) 



^The Landau gauge condition in momentum space, q^Df^^{q) = places a constraint on the 
Lorentz components of the propagator so that, for non-zero momentum, there are — 1 degrees 
of freedom [14|. 
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and 



m = ^T.T.DM- (2.11) 



Nfi is the number of space-time dimensions and the available momentum values, q, are given 
by 

g. = ^, n,e{--,-\. (2.12) 

The range of q is determined by the fact that our lattices have an even number of points in 
each direction and that we use periodic boundary conditions. In the continuum, the scalar 
propagator is related to the full propagator by 

D;'M = (V - '^-^W'D{q^) (2.13) 

in Landau gauge. 

Landau gauge is a smooth gauge that preserves the Lorentz invariance of the theory, so 
it is a popular choice. We work in Landau gauge for ease of comparison with other studies, 
and because it is the simplest covariant gauge to implement on the lattice. All configurations 
were gauge fixed by maximizing an O{o^) improved Landau gauge fixing functional using 
Conjugate Gradient Fourier Acceleration [TH] as described in Ref. [jl6|. 



III. TREE-LEVEL CORRECTION 

One thing that is known about the gluon propagator is its perturbative, asymptotic 
behavior. In the spirit of improvement, we can use this knowledge to augment our lattice 
results and make better contact with the continuum. In the continuum, as — oo, the 
propagator has the form 

m = \ (3.1) 

up to logarithmic corrections. A well known artifact of the lattice is that for a free massless 
boson with an unimproved action the lattice propagator has the form 



It has been argued, in Ref. |]T^ and elsewhere, that the correct momentum variable to 



use when examining the gluon propagator on the lattice, with the Wilson action, is not 
Eq. (|12D, butg 



^Many authors have q and q defined the other way around, but in this context our terminology is 
more natural. 
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_ 2 . %a 

q^ = -sm—. (3.3) 

It has been observed that this choice ensures that the propagator takes its asymptotic form 
at large lattice momenta p7| , p^ . 

The improved action Eq. ( |2.1| ) together with the gluon field defined in Eq. ( p.6| ) has the 
C(a^) improved tree- level behavior []^,||] 

= (3.4) 



and we will use Eqs. ( p.2|) and ( p.4|) to obtain the correct momentum variable for each 
action. To emphasize the nonperturbative aspects of the propagator, we divide it by its 
perturbative, 1/g^ result. Hence, all figures are plotted against q^D{q^), which is expected 
to approach a constant up to logarithmic corrections as — > oo. We will see that this also 
makes for a stringent test of the ultraviolet behavior of the propagator. We will work with 
the momentum variables defined as 



g;r-^sin^, (3.5) 



and 



for the Wilson and improved actions respectively. A similar momentum variable was used 
in the study of the gluon propagator in Ref. . 



In the language of continuum physics 

where n(p^) is the scalar vacuum polarization. In the asymptotic region, 1/[1 + n(j9^)] 1 
up to logarithmic corrections. We argue that it is the lattice version of D{p'^) / D^^'^^ijp) 
that will most rapidly approach its continuum form as the lattice spacing is reduced and 
we will later graphically demonstrate this. The essential point is that at large momenta the 
lattice gluon propagator will experience asymptotic freedom just as in the continuum, i.e., 
the ultraviolet propagator will approach its tree-level form. Thus on the lattice we expect 
to find D{p'^) / D^^^'^{p'^) — i> 1 for large even though the ultraviolet lattice artifacts in both 
DijP') and D^^^°{p^) may themselves be large. We will refer to this procedure for minimizing 
ultraviolet lattice artifacts as tree-level correction. This philosophy is similar to that applied 
in recent studies of the quark propagator pO[. In figures where q^D{q^) is plotted vs. g, the 
"g" in q'^D{q^) (plotted on the y-axis) is always the same as the q that is used on the x-axis, 
where q = q, q^ or q^ as described in the text. 

The bare, dimensionless lattice gluon propagator D{qa) is related to the renormalized 
continuum propagator Dn^q; /i) by 

a^D{qa) = Z^i^fi, a)DR{q; fi), (3.8) 
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Dimensions 


/3 


a (fm) 


Volmne (fm'^) 


Configurations 


Iw 


16^ X 32 


5.70 


0.179 


2.87^ X 5.73 


100 


li 


16^ X 32 


4.38 


0.166 


2.64^ X 5.28 


100 


2 


10^ X 20 


3.92 


0.353 


3.53^ X 7.06 


100 


3 


8^ X 16 


3.75 


0.413 


3.30^ X 6.60 


100 


4 


16=^ X 32 


3.92 


0.353 


5.65^ X 11.30 


100 


5 


12^ X 24 


4.10 


0.270 


3.24^ X 6.48 


100 


6 


32^ X 64 


6.00 


0.099 


3.18^ X 6.34 


75 



TABLE I. Details of the lattices used to calculate the gluon propagator. Lattices Iw and li 
have the same dimensions and approximately the same lattice spacing, but were generated with 
the Wilson and improved actions respectively. Lattice 6 was generated with the Wilson action. 



for momenta, q, sufficiently small compared to the cutoff, n/a. D^^q] jj) is independent of a 
for sufficiently fine lattices; i.e. in the scaling regime. The renormalization constant Z^{ix, a) 
is determined by imposing a renormalization condition at some chosen renormalization scale 
/i, e.g., 

DR{q%.^^. = \. (3.9) 

The renormalized gluon propagator can be computed both nonperturbatively on the lattice 
and perturbatively in the continuum for choices of the renormalization point in the ultravi- 
olet. It can then be related to the propagator in other continuum renormalization schemes 
such as MS. 



IV. RESULTS 



A. Analysis Overview 



The gluon propagator has been calculated on seven different lattices, the details of which 
are listed in Table I. Note that the ffist two are labeled "Iw" and "li" . These have the same 
number of lattice points at almost the same spacing (hence approximately the same physical 
volume), but Iw was generated with the standard, Wilson gauge action, while li used the 
C(a^) improved action ( p.l|) . Lattice 6 was generated with the Wilson action and used to 
study the gluon propagator in Ref . |]r7| . A value for the tadpole factor has been obtained for 
(3 = 6.0 of uq = 0.878 and this has been used to normalize the propagator with respect to 
the other lattices. It will be used here for comparison purposes as it is finer than the other 
lattices. Configurations on lattices 2-5 were generated with the O(a^) improved action. All 
of the propagators are plotted in physical units, where the scale has been determined by the 
static quark potential with a string tension of ^/a = 440 MeV. Details of this calculation 
may be found in Ref. . 

Data points that come from momenta lying entirely along a spatial Cartesian direction 
are indicated with a square while points from momenta entirely in the temporal direction are 
marked with a triangle. As the time direction is longer than the spatial ones any difference 



6 



between squares and triangles may indicate that the propagator is affected by the finite 
volume of the lattice. Data points from momenta entirely on the four-diagonal are marked 
with a diamond. Systematic separation of data points taken on the diagonal from those in 
other directions indicates violation of rotational symmetry. 

In the continuum, the scalar function is rotationally invariant. Although the hypercubic 
lattice breaks 0(4) invariance, it does preserve the subgroup of discrete rotations Z(4). In 
our case, this symmetry is reduced to Z(3) as one dimension will be twice as long as the 
other three in each of the cases studied. We exploit this discrete rotational symmetry to 
improve statistics through Z(3) averaging. This is best explained through a simple example. 
Consider the propagator at momentum q = (3, 2, 1,4) (say). Z(3) symmetry means that 

D(3, 2, 1, 4) = D{2, 3, 1, 4) = D{2, 1, 3, 4) = D{1, 2, 3, 4) = D(l, 3, 2, 4) = D(3, 1, 2, 4) (4.1) 

so we calculate the propagator for each of these values of momentum, and then average the 
results. 



B. Tree-Level Correction and Rotational Symmetry 

The "raw" gluon propagator from lattices Iw and li is shown in Figures |l] and |] re- 
spectively. Both of these have been plotted as functions of q, Eq. ( p.l2|) , for all available 



momenta, and both show severe ultraviolet noise. We may take some comfort from the 
observation that the signal degradation is not as bad in the improved case where the finite 
spacing errors do not exceed the infrared peak and the UV tail is generally flatter. However, 
neither result looks at all satisfactory at large momenta. No data cuts or tree-level correction 
have yet been used. 

The most obvious way to deal with this noise is to apply an ultraviolet cut, considering 
only momenta out to half of the Brillouin zone. For each of the four Cartesian directions, 

q < f • (4.2) 

2a 

We refer to this as the "half-cut" and in Fig. ^ and Fig. |^ we see that this removes the worst 
of the artifacts. The two propagators, show plausible asymptotic behaviors, but there are 
still clear signs of lattice artifacts and we have lost a lot of data in the ultraviolet. While 
neither of these shortcomings is a significant problem for studies of the infrared, we will 
show that something as crude as the half-cut is not necessary and we can do much better 
at minimizing lattice artifacts. 

We have already argued the case for applying a tree-level correction through the use 
of the alternative momentum variables derived from the tree-level behavior of the actions. 
The effect of doing this is seen in Fig. |^ and Fig. ^, where the Wilson propagator has been 
plotted as a function of q^ {q^D{q^) vs. q^) and the improved propagator as a function 
of q^ {q^D{q^) vs. q^) for all momenta of the Brillouin zone. Comparing these to Figs. |l| 
and I], we see an excellent restoration of rotational symmetry all the way to the edge of the 
Brillouin zone. This is especially true of the improved action case in Fig. |^. The propagators 
also appear to be approaching their asymptotic, perturbative values. Later, momentum cuts 
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FIG. 1. Uncorrected gluon propagator from lattice Iw (/? = 5.70,16^ x 32, Wilson action), 
plotted as a function of q. The dramatic "fanning" is caused by finite spacing errors which quickly 
destroy the signal at large momenta. 
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FIG. 3. Uncorrected gluon propagator from lattice Iw (/3 = 5.70,16^ x 32, Wilson action), 
plotted as a function of q with the momentum "half-cut" appHed. 
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FIG. 4. Uncorrected gluon propagator from lattice li (/3 = 4.38,16^ x 32, improved action), 
plotted as a function of q with the momentum "half-cut" applied. The improved propagator has 
different normalization to the Wilson case due to a difference in the Z3 renormalization constant. 
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FIG. 5. Uncut gluon propagator from lattice Iw (/3 = 5.70,16^ x 32, Wilson action), plotted 
as a function of for all momenta. The tree-level correction has greatly reduced discretization 
errors from those seen in Fig. |l[ 
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FIG. 6. Uncut gluon propagator from lattice li (/? = 4.38, 16^ x 32, improved action), plotted 
as a function of for all momenta. The combination of improved action and tree-level correction 
has produced a remarkably clean signal over the entire range of accessible momenta. This figure 
should be compared with Fig. ||, and with Fig. ^ for the Wilson action at a similar lattice spacing. 
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will be applied to the data to further eliminate lattice artifacts, but for the moment it is 
interesting to keep all data, as they provide insight into the behavior of lattice simulations. 

Both Figs. H and |^ are consistent with the study of Ref . [l^ , but the discrepancy between 
diagonal and Cartesian points in Fig. |^ is a clear sign of rotational symmetry breaking in 
the unimproved case. With the Wilson action, the quality of the data is suffering from 
the coarseness of the lattice. As we might hope, the improved propagator in Fig. |^ shows 
excellent agreement between diagonal and Cartesian points, and the data is generally less 
spread. The propagator from the improved action has better rotational symmetry at the 
same lattice spacing. Less easy to understand is the slight suppression of the temporal 
points (triangles) in the Wilson case. Fig. ^. The time axis of this lattice (as with all the 
lattices considered here) is twice as long as the other three axes, so different values for the 
points along the long axis would normally be interpreted as a finite volume effect, yet there 
is no sign of it in the improved case (which has approximately the same physical volume). 
There is a difference between the improved and unimproved cases in the amplitudes of the 
propagators, but this is accounted for by renormalization and will be discussed below. 

Out of curiosity the gluon propagator from lattice li has also been examined as a function 
of g^, which we have already argued to be inappropriate. Not surprisingly, this leads to a 
"propagator" that suffers badly from lattice artifacts. We have not included a figure here, 
but the resulting propagator droops strongly in the ultraviolet. This is clearly a poor choice 
of momentum variable for this action as expected on the basis of our tree-level correction. 
For best results at finite lattice spacing, the correct momentum variable is determined by 
the appropriate tree-level behavior, which in turn is defined by the choice of action and 
gluon field definition. For the rest of this report it shall be implicit that when discussing 
quantities from the Wilson action, is used, and is used with the improved action. 



C. Lattice Spacing Dependence 

At this point it is interesting to explore the effect of making the lattice coarser. Fig- 
ures 1^, I and ^ show the uncut, tree-level corrected propagator on progressively coarser 
lattices (a = 0.27, 0.35 and 0.41 fm respectively). Consider the most extreme case, shown in 
Fig. ^ This very coarse lattice has spacing a = 0.41 fm, which is more than twice as coarse 
as the previous lattices. Any sign of a perturbative tail has been lost, as the UV cutoff has 
been lowered, but the infrared behavior remains. There is no sign of any qualitative change, 
which appears to indicate that even on such a coarse lattice we are not losing information 
vital to the infrared physics of the gluon propagator. 

This gives us great confidence in the use of improved actions on coarse lattices for the 
probing of nonperturbative physics. This is the motivation for creating lattice 4 at a = 0.35 
on a very large volume. Fig [l^, which shows the results from this large lattice, shows no 
signs of significant finite volume artifacts when compared with Fig. |^ which has the same 
lattice spacing, but a smaller volume. 
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FIG. 7. Gluon propagator from lattice 5 at /3 = 4.10, which has spacing a ~ 0.27 fm on 12^x24. 
This has the same physical volume as lattice 3 of Fig. ^. The propagator is shown for all momenta 
(no data cuts) after tree-level correction. 




1 - 



1 2 3 4 5 6 

q (GeV) 

FIG. 8. Gluon propagator from lattice 2, the smaller lattice at /3 = 3.92 which has spacing 
a ~ 0.35 fm on a 10^ x 20 lattice. Finite volume errors are just detectable as indicated by momenta 
along the time axis (filled triangles) falling below the rest of the data. Tree-level correction has 
been used, but no data cuts have been applied. 
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FIG. 9. Gluon propagator from lattice 3 at /? = 3.75, which has spacing a ~ 0.41 fm on 8^ x 16. 
The propagator is shown for aU momenta (no data cuts) after tree-level correction. This propagator 
is consistent with that obtained on much finer lattices. 




1 2 3 4 5 6 

q (GeV) 

FIG. 10. Gluon propagator from lattice 4, the larger lattice at (3 = 3.92, which has spacing 
a ~ 0.35 fm on a 16'^ x 32 lattice providing the largest physical volume of any in this study. 
Tree-level correction has been used, but no data cuts have been applied. 
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D. Data Cuts 



Having identified possible lattice artifacts, cuts may be applied to clean up the data, 
making it easier to draw conclusions about continuum physics. Data at large momenta will 
of course be most susceptible to finite lattice spacing errors. We choose to prefer data from 
momentum points near the four- diagonal, as this evenly samples all Cartesian directions, 
i.e., for a given momentum squared (g^) it has the smallest values of each of the Cartesian 
components q^. This should minimize finite lattice spacing artifacts. 

We calculate the distance Ag of a momentum vector q from the diagonal using 

Ag = |g| sin6'g, (4.3) 

where the angle 6q is given by 

cos^, = (4.4) 

and n = |(1, 1, 1, 1) is the unit vector along the diagonal. In this way we ignore data points 
that are potentially most affected by hypercubic artifacts. We call this cut the cylinder 
cut [jl^. From this point on, we exclude points greater than two spatial momentum unitsQ 



from the four-diagonal. Furthermore, the point at zero four-momentum has been cut from 
all the following plots of q^D{q). On any finite lattice, -D(O) must be finite, hence q^D{q) = 
for q = 0. This point is therefore trivial when plotting q'^D{q). When the scalar function, 
D{q), itself is considered we can make a study of -D(O) by considering it on lattices of 
differing volumes and then making an infinite volume extrapolation. We will perform this 
extrapolation below. 



E. Action Dependence 

Once again we compare the gluon propagator generated with the Wilson action to that 
generated with the improved action after tree-level correction, this time applying the cylin- 
der cut. To make the comparison in Fig. 0, we note that there is of course a small difference 
in normalization. This is the difference in the renormalization between the Wilson and 
improved propagators. As the relative renormalization is independent, the unimproved 
propagator has been multiplied by a relative renormalization of 1.09 to make direct com- 
parison possible. This factor is deduced by adjusting the vertical scales of the two data sets 
until they agreed. Apart from the superior performance of the improved propagator, which 
has already been discussed, the two actions produce the same result. 

We push our results further by comparing the improved f3 = 4.38 propagator with that 
from lattice 6 (Wilson action), which is finer (a = 0.1 fm), has more points (32'^ x 64) and is a 
little larger. Both data sets are cylinder cut, and each is tree-level corrected according to its 



spatial momentum unit is 2tt /aLg where Ls is the number of lattice sites in the spatial directions 
{Lg = = Ly = Lz). 
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action. The relative renormalization has been determined to be improved )/Z3( Wilson) 



1.08. It can be seen from Fig. ^ that not only are the two propagators consistent, but 



that the ultraviolet performance of lattice li is remarkable. The propagator from Ref. |T7 



had the momentum half-cut applied, whereas our improved propagator with lattice spacing 
a = 0.17 fm is shown for the entire Brillouin zone. We have calculated the propagator over 
the same range of momenta as Ref. Jl^, despite using a much coarser lattice. 



F. Scaling Analysis 



Next, we consider the propagator on the coarser lattices. Fig. |13| shows the propagator 
from lattices li, 2, 3 and 5. Examining Figures |TT| and |TB| we see that the Wilson f3 = 5.7 
and improved (3 = 4.10 and P = 4.38 results all agree well, which suggests that these are 
"fine enough" lattices. We see that the /3 = 3.75 and P = 3.92 propagators do not quite hne 
up with the others, but instead the UV tail rises slightly as the lattice becomes coarser. This 
is an indication of a loss of scaling. The lattices at /5 = 3.92 and (3 = 3.75 having a = 0.35 
and 0.41 fm respectively are too coarse for the tree-level correction to completely correct 
the entire Brillouin zone, which is not surprising. We have placed extraordinary demands 
on our simulations by examining them near the cutoff. The conclusion is that such coarse 
lattices should be half-cut. Nevertheless, the propagators all agree in the infrared. Now that 
we have an understanding of the dependence of lattice propagator on the lattice spacing, we 
can study the effect of the finite volume. 

G. Volume Dependence 



Results from lattices 2 and 4 have already been reported in Ref. |]I0[ and are presented 
again here for completeness and ease of comparison. They have same lattice spacing, but 
different numbers of lattice points, and hence different physical volumes. The gluon prop- 



agator has been calculated on each lattice, and the results compared in Fig. |TJ. The two 
propagators are consistent in this figure, despite the fact that one lattice has sides 60% 
longer in all four directions. This shows that finite volume effects are small compared to the 
statistical errors. The turn over seen in the gluon propagator in lattice studies is certainly 
not a finite volume effect. Note that 5.65^ x 11.30 fm^ is a very large volume by the stan- 
dards of present day lattice studies, and gives us an unprecedented look at the behavior of 
QCD in the deep infrared. 

Fig. shows the cylinder-cut data for the scalar function D{q^) for each of the improved 
lattices. This plot provides a dramatic demonstration of lattice artifacts. In this way of 
plotting our results, the five lattices appear in very good agreement in the ultraviolet and 
through intermediate momenta. When plotted in this way, we can see that below ~ 500 MeV 
the propagators do begin to differ due to finite volume effects. As the volume increases, the 
low momenta data points drop, until we can see the infrared flatten off. The grouping of 
points around 400 MeV suggest that we have, for the two largest lattices, results indicative 
of the infinite volume limit. At ~ 250 MeV, the results for the two largest lattices (both 
/3 = 3.92) are consistent, and in particular the fact that the small difference between them 
is produced by such a large difference in volume gives us confidence in the results. For 
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comparison, the tree-level, perturbative expression D{q'^) = is also shown, suitably 
normalized. 

It is interesting to note that the disagreement in the propagators above 1 GeV or so 
revealed in Fig. |l^ is hidden by the scale of the vertical axis in Fig. |15|. Multiplication of 
the propagator by is required to amplify this region and critically examine the extent to 
which lattice spacing artifacts are removed by O(a^) improvement terms. A failure to do 
this could lead to incorrect conclusions being drawn on the effectiveness of improvement in 
the gluon propagator. Thus it is always best to plot q^D{q^) versus g^, when the hypercubic 
artifacts are of interest. 



H. Asymptotic Behavior 

For further comparison with perturbation theory, we have chosen to show the gluon 
propagator from 1.5 to 5.5 GeV, in Fig. |TB|. In this window, the transition from perturbative 
to nonperturbative physics can be clearly seen. As well as the lattice gluon propagator and 
the tree-level, continuum propagator, we show a perturbative, three-loop calculation 



We used parameters obtained from Ref. [g], where at the renormalization point, /i = 5.48 
GeV, the strong coupling constant was found to be = 0.255. That was a quenched 
calculation, so this number should not be compared directly with experiment. The data 
agree very well with three-loop perturbation theory down to g ~ 2.5 GeV. Below 2 GeV we 
see that three-loop perturbation theory begins to fail. 



I. Propagator at Zero Four-Momentum 

Values for the gluon propagator at zero four-momentum are shown in Table II for each 
of the lattices created in this investigation. Statistical errors are given in parentheses. The 
renormalization condition of Eq. ( |3.9D is enforced at the renormalization point = 4.0 
GeV, which sets the scale for D{q^). We see that as the volume of the lattice increases, 
Z^(0) becomes smaller. In Fig. |l^ we plot the infrared behavior of the renormalized gluon 
propagator for five lattices and we include the points calculated at zero momentum in this 
plot. We see that the the infrared behavior is quite smooth and reasonably consistent for 
our two largest volume lattices {(5 = 3.92, small and large). Fig. |T^ illustrates the data with 
a linear fit in the inverse volume according to 

DiO) = + Do^iO) (4.5) 

We find a reasonable fit with parameter values c = 245(22) fm^ GeV~^ and Dao{0) = 7.95(13) 
GeV~^, where Doo(O) is the infinite volume limit of the zero-momentum gluon propagator. 
Fig. ITBI strongly supports the hypothesis that the gluon propagator is finite in the infrared. 



It is also clear that the results of our largest physical volume lattice are very close to the 
infinite volume limit. 

Note that a complete systematic extrapolation to the infinite volume limit remains to be 
carried out in the future. Ideally, one performs a number of calculations at fixed volume and 
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0.0 ' ^ ^ ^ ^ ^ ' 

1 2 3 4 5 6 
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FIG. 11. Comparison of the gluon propagator from lattices Iw at /? = 5.70 and li at /? = 4.38. 
Data has been cylinder cut and tree-level correction has been applied. We have determined 
(improved) /Z3 (Wilson) = 1.09 by matching the vertical scales of the data. 



Lattice 


Dimensions 


P 


D(0) 


D(0) (GeV-2) 


Volume (fm^) 


li 


16^ X 32 


4.38 


32.0 (8) 


10.4 (2) 


97.2 


Iw 


16^ X 32 


5.70 


24.0 (5) 


10.0 (2) 


135 


5 


12^ X 24 


4.10 


10.6 (3) 


9.0 (2) 


220 


3 


8^ X 16 


3.75 


4.3 (1) 


8.9 (2) 


237 


2 


10^ X 20 


3.92 


5.7 (1) 


8.6 (2) 


300 


4 


16^ X 32 


3.92 


5.4 (1) 


8.2 (2) 


2038 



TABLE IL The value of gluon propagator at zero four-momentum for each of the lattices 
created in this investigation, in order of increasing volume. The raw (dimensionless) and physi- 
cal values are given. In obtaining the physical values we have set the renormalization condition 
D(/x^) = l//i^ at /X = 4.0 GeV. An estimate of the uncertainty in the last figure is given in 
parentheses. 
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FIG. 12. Comparison of the gluon propagator from the finest improved lattice (lattice li, 
P = 4.38) and the finest Wilson lattice (lattice 6, f3 = 6.0). Data has been cylinder cut and the 
appropriate trcc-lcvcl corrections have been applied. The data from lattice 6 is half-cut whereas 
lattice li displays the full Brillouin zone. We have determined Z3 (improved)/Z3 (Wilson) = 1.08 
by matching the vertical scales of the data. 
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FIG. 13. Comparison of the gluon propagator from lattices li {(3 = 4.38), 2 (/3 = 3.92, small), 
3 (/3 = 3.75), and 5 (/3 = 4.10), which have a variety of lattice spacings. Data has been cylinder 
cut and tree- level correction has been applied. Data from the two finest improved lattices (0.17 
and 0.27 fm) are consistent. A clear violation of scaling is seen in the coarsest two lattices (0.35 
and 0.41 fm), where the spacing is too coarse for tree-level correction to completely restore the full 
Brillouin zone behavior. 




3 

q (GeV) 

FIG. 14. Comparison of the gluon propagator from lattices 2 and 4, which have the same lattice 
spacing (a = 0.35 fm), but different lattice volumes. Notice that one lies directly over the other, 
despite having very different volumes. Data has been cylinder cut and tree-level correction has 
been applied. 
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FIG. 15. Comparison of the gluon propagator generated with an improved action on five differ- 
ent lattices. We find good agreement down to g ~ 500 MeV. At the lowest accessible momenta the 
data points drop monotonically with increasing volume, but the lowest point (on the largest lattice) 
shows signs of having converged to its infinite volume value. For comparison with perturbation 
theory, a plot of the continuum, tree-level gluon propagator (i.e., Xjq^ appropriately scaled) has 
been included. 



various lattice spacings and then performs a continuum limit extrapolation for that fixed 
volume. This continuum limit extrapolation would be done for each of a variety of lattice 
volumes and then finally an infinite volume extrapolation performed on those results. This 
procedure corresponds to the axiomatic field theory prescription of taking the continuum 
limit before the infinite volume limit. Given this caution, the finite precision of this study, 
and the fact that the linear ansatz above may be incorrect, we can not completely exclude 
the possibility that the deep infrared (i.e., below ~ 350 MeV) behavior of D{q^) may very 
slowly decrease toward zero as the infinite volume limit is taken. 

It is interesting to compare our results with a recent calculation of the gluon propagator 
in Laplacian gauge p3|, which is expected to be free of gauge ambiguity. In that gauge, the 
propagator takes its perturbative. Landau-gauge value in the asymptotic region and is also 
infrared finite. The Laplacian gauge propagator is seen to have a behavior similar to that 
seen here. 



V. CONCLUSIONS 

The gluon propagator has been calculated on a set of lattices with an O(a^) mean- field 
improved action, in 0{a?) mean-field improved Landau gauge. Tree-level correction has been 
shown to reduce rotational symmetry breaking and dramatically improve the ultraviolet 
behavior of the propagator. 
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FIG. 16. Comparison of the lattice gluon propagator with that obtained from perturbation 
theory, in the ultraviolet to intermediate regime. The continuum expressions are tree-level (i.e., 
1/g^ appropriately scaled) and the three- loop expression used in Ref . [2] . 

For /3 > 4.10 (a < 0.27 fm) the tree- level corrected improved propagator displays scaling 
over the entire Brillouin zone. Kt (3 = 4.38 (a = 0.166 fm), the gluon propagator has excellent 
behavior for the entire range of available momenta in the Brillouin zone, reproducing the 
anticipated UV behavior of perturbation theory to three-loops. 

The infrared behavior of the gluon propagator is robust even with a lattice spacing of 
0.41 fm. Calculation on a lattice with a large volume indicates that finite volume effects 
are small. In particular, the turn over observed in previous studies of the Landau gauge 
gluon propagator is not a finite volume artifact. We conclude that the propagator is almost 
certainly infrared finite, in agreement with earlier studies. A significant volume dependence 
is revealed only at the smallest non-trivial momenta. An extrapolation of D{Q) via a linear 
ansatz inversely proportional to the physical lattice volume provides a reasonable fit. More- 
over, results from our largest volume lattice reside very close to the infinite volume limit. 
We have probed the approach to the infinite volume limit by first determining a range of j3 
in which the propagator scales for q < 0.7 GeV on similar finite volumes. Physically large 
volumes are accessed by decreasing (3 within the scaling range on large lattices. A more 
complete study of the infinite volume limit should be undertaken in the near future. 

The tree-level corrected results from our 13 = 3.92 (a = 0.353 fm) 16'^ x 32 lattice with 
a physical volume of 5.65^ x 11.30 = 2038 fm^ may be regarded as an excellent estimate 
of the infinite volume, continuum limit Landau-gauge gluon propagator for q < 0.7 GeV. 
The tree-level corrected results from our /3 = 4.38, (a = 0.166 fm) results presented here 
are an excellent estimate of the infinite volume, continuum limit of the Landau-gauge gluon 
propagator for q > 0.7 GeV. We have seen that these two sets of data smoothly match in 
the intermediate regime {q ~ 0.7 GeV) and are entirely consistent with each other in this 
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FIG. 17. The renormalized gluon propagator is shown in the infrared region, including the 
points at zero four-momentum, from five lattices. 
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region. The possible effects of lattice Gribov copies remains a very interesting question and 
we plan to extend this study to Laplacian gauge and other related gauge-fixing schemes in 
the near future. 
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FIG. 18. Values for the gluon propagator at zero four-momentum, -D(O), plotted as a function 
of the inverse lattice volume. The solid line represents a linear fit to the lattice results. The fit 
indicates the largest volume results are very close to the infinite volume limit and -D(O) = 7.95(13) 
GeV~^ in the infinite volume limit. 
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